Mori. Not. R. Astron. Soc. 000, 1-8 (0000) Printed 27 February 2013 (MN I^TeX style file v2.2) 

From light to mass: accessing the initial and present-day 
Galactic globular cluster mass functions 



(N 

o 



C. Bonatto 1 and E. Bica 1 

1 Departamento de Astronomia, Universidade Federal do Rio Grande do Sul, Av. Bento Goncalves 9500 
Porto Alegre 91501-970, RS, Brazil 



27 February 2013 



o 

(N 

< 

o 

6 



> 

in 

en 

o 

(N 



>i-( ' 1 INTRODUCTION 



3 



ABSTRACT 

The initial and present-day mass functions (ICMF and PDMF, respectively) of the 
Galactic globular clusters (GCs) are constructed based on their observed luminosities, 
the stellar evolution and dynamical mass-loss processes, and the mass-to-light ratio 
(MLR). Under these conditions, a Schechter-like ICMF is evolved for approximately a 
Hubble time and converted into the luminosity function (LF), which requires finding 
the values of 5 free parameters: the mean GC age (t gc ), the dissolution timescale of a 
10 5 Mq cluster (rg), the exponential truncation mass (M c ) and 2 MLR parametris- 
ing constants. This is achieved by minimising the residuals between the evolved and 
observed LFs, with the minimum residuals and realistic parameters obtained with 
MLRs that increase with luminosity (or mass) . The optimum PMDFs indicate a total 
stellar mass of ~ 4 x 10 7 M Q still bound to GCs, representing ~ 15% of the mass in 
clusters at the beginning of the gas-free evolution. The corresponding ICMFs resem- 
ble the scale-free MFs of young clusters and molecular clouds observed in the local 
Universe, while the PDMFs follow closely a lognormal distribution with a turnover 
at Mto ~7x 10 4 M Q . For most of the GC mass range, we find an MLR lower than 
usually adopted, which explains the somewhat low Mto- Our results confirm that the 
MLR increases with cluster mass (or luminosity), and suggest that GCs and young 
clusters share a common origin in terms of physical processes related to formation. 

Key words: (Galaxy:) globular clusters: general 



Globular clusters (GCs) are the living fossils of an epoch 
dating back to the time of Milky Way's formation and thus, 
they probably share the same physical conditions and pro- 
cesses that originated our Galaxy. Since formation, GCs are 
continuously affected by several sources of mass loss but, be- 
cause of a large initial mass and/or an orbit that takes them 
most of the time away from the Galactic centre and/or disk, 
they are, in general, characterised by a long longevity. This 
suggests that some primordial information may have been 
preserved, if not at the individual scale, at least with re- 
spect to the collective properties shared by the surviving 
GC population. 

In this context, an interesting issue is the origin of the 
Milky Way's (old) GC population and its relation to the 
young massive clusters that form in starbursts and galaxy 
mergers of the local Universe. If both old and young cluster 
populations share a similar, starburst origin, the mass dis- 
tribution of their members should also be characterised by 
similar mass functions (MFs), something that is ruled out by 
observations in several local galaxies. For instance, when ex- 



pressed in terms of the number of GCs per unit logarithmic 
mass- interval (dN/dlog M), the MFs of young clusters in 
several galaxies follow monotonically a power law of slope 
» -1, while the present-day GC MF (PDMF) presents a 
turnover for masses lower than Mto ~ 10° Mq (McLaughlin 
& Fall 2008 and references therein). The low-mass drop in 
the lognormal PDMF has been interpreted as mainly due 
to the large-scale dissolution of low-mass GCs after several 
Gyrs of tidal interactions with the Galaxy and internal dy- 
namical processes (e.g. McLaughlin & Fall 2008; Kruijssen 
2008; Kruijssen & Portegies Zwart 2009). 

Consequently, since mass loss continuously changes the 
MF shape - especially in the low-mass range, a direct com- 
parison of the several Gyr-old PDMF with the MF of young 
clusters should not be used to draw conclusions on a phys- 
ical relation between both populations. On the other hand, 
this raises another, but related question of whether or not 
the initial MFs (ICMFs) of globular and young clusters are 
similar, since theoretical evidence suggest that cluster for- 
mation leaves distinct imprints in the ICMFs (Elmegreen 
& Falgarone 1996). Thus, if both cluster populations have 
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the same origin, this should be reflected in the respective 
ICMFs. 

Given the above, the reconstruction of Milky Way's GC 
ICMF requires having access not only to the (as complete 
as possible) PDMF, but also to the mass-loss physics. As for 
the latter, recent numerical simulations and semi-analytical 
models have produced robust relations for the mass evolu- 
tion of star clusters located in different environments and 
losing mass through stellar evolution and tidal interactions 
with Galactic substructures and giant molecular clouds (e.g. 
Lamers et al. 2005; Lamers, Baumgardt & Gieles 2010). 

Over the years, the GCMF time evolution has been in- 
vestigated with a wide variety of approaches. For instance, 
Vesperini (1998) takes the effects of stellar evolution, two- 
body relaxation, disc shocking, dynamical friction, and the 
location in the Galaxy into account to follow the evolution 
of individual GCs. The main findings are: the large-scale dis- 
ruption of GCs in the central parts of the Galaxy leads to a 
flattening in their spatial distribution; if the initial GCMF 
shape is lognormal, it is preserved in most cases; and in the 
case of an initial power-law GCMF, evolutionary processes 
tend to change it into lognormal. Baumgardt (1998) consid- 
ers GC systems formed with a power-law MF and follows 
the orbits of individual clusters. Taking several GC destruc- 
tion mechanisms into account, he finds luminosity distribu- 
tions that depend on Galactocentric distance. Vesperini et 
al. (2003) study the dynamical evolution of M 87 GCs with 
numerical simulations that cover a range of initial conditions 
for the GCMF and for the spatial and velocity distribution 
of the GCs. The simulations include the effects of two-body 
relaxation, dynamical friction, and stellar evolution-related 
mass loss. They confirm that an initial power-law GCMF 
can be changed into bell-shaped through dynamical pro- 
cesses, and show that the inner flattening observed in the 
spatial distribution of M 87 GCs can result from the dynam- 
ical evolution on an initially steep density profile. Parmen- 
tier & Gilmore (2007) show that an approximately Gaussian 
MF is naturally generated from a power-law mass distribu- 
tion (but with a lower-mass limit) of protoglobular clouds 
by expulsion from the protocluster of star-forming gas due 
to supernova activity. Shin, Kim & Takahashi (2008) tar- 
gets the MF and radial distribution evolution of the Galac- 
tic GCs by means of a Fokker-Planck model that considers 
dynamical friction, disc/bulge shocking and eccentric GC 
orbits. Their best-fit models predict the initial GC mass as 
Mqc ~ (1-5-1.8) x 10 s M Q . 

At some point, most of the previous works compare 
models with the observed GC PDMF. However, especially 
because of the large distances involved, a direct determina- 
tion of the stellar mass for most GCs is not available, which 
means that the observational PDMF cannot be taken as an 
ideal starting point for this analysis. Instead, the most direct 
and robust observable available for the present-day Milky 
Way's GC population as a whole is the luminosity function 
(LF), which is usually expressed in terms of the absolute 
V magnitude. Obviously, the use of luminosity requires a 
mass-to-light relation (MLR) to convert LFs into MFs. 

Recent observational and theoretical findings have con- 
vincingly shown that the MLR is not constant with respect 
to cluster mass. Instead, the dynamical evolution of clusters 
- and the preferential loss of low-mass stars - also affects the 
MLR. Specifically, for a given age and metallicity, the MLR 



is expected to increase with cluster mass (or luminosity; e.g. 
Kruijssen 2008; Kruijssen & Portegies Zwart 2009). Inter- 
estingly, the MLR dependence on luminosity (or mass) is 
already suggested by the dynamical cluster mass and MLR 
estimates of Pryor & Meylan (1993). 

Our goal with the present paper is to reconstruct both 
the ICMF and PDMF of the Galactic GCs, starting with 
their observed LF and taking the mass loss into account. In 
this process, different analytical forms of the MLR will be 
tested to search for constraints. 

The present paper is organised as follows. In Sect. 2 
we build the observed and theoretical LFs of the Galactic 
GCs, convert the latter into the PDMF, and discuss the 
adopted MLRs. In Sect. 3 we briefly describe the parameter- 
search method. In Sect. 4 we discuss the results. Concluding 
remarks are given in Sect. 5. 



2 THE THEORETICAL PRESENT-DAY LF 

The first step in our analysis is to build the observed (or 
present-day) LF (<j)(M v ) = dN/dM v , i.e., number of GCs 
per unit absolute magnitude interval) of the Galactic GCs. 
The data were taken from the 2010 revision of the catalogue 
compiled by Harris (1996), which contains 156 GCs with 
bona-fide absolute magnitudes in the V band (My)- By def- 
inition, the integral of (j>(Mv), over the full Mv range is 
the input number of GCs. As a caveat, we note that several 
works indicate that not all Galactic GCs share a common 
history (e.g. Bica et al. 2006 and references therein), with 
some having been accreted from dwarf galaxies (e.g. Hilker 
& Richtler 2000; Parmentier et al. 2000) and others associ- 
ated with the Galactic thick disc (e.g. Zinn 1985). 

Following previous work (e.g. Bonatto & Bica 2011), 
uncertainties in My are incorporated in the LF by assum- 
ing a normal (Gaussian) distribution of errors. Specifically, 
if the magnitude and standard deviation of a GC (usually, 
the mean value over a series of independent measurements) 
are given by Myiu, the probability of finding it at a specific 

! / M v - My \ 2 

value M v is given by P{M V ) = -?=; e 2 V * ) . Next, 
we build a grid of bins covering the full Mv range in Harris 
(1996) and, for each bin, we sum the Mv density (or prob- 
ability) for all GCs (i.e., the difference of the error function 
computed at the bin borders). However, since uncertainties 
are not given in Harris (1996), for simplicity we adopt rel- 
ative errors that depend on the apparent magnitude (mv) 
as e = 0.1e ( - mv/27) . Thus, relatively bright (mv < 10) GCs 
have e < 0.15, while fainter ones (mv > 20) have e > 0.2. 
Finally, we assume that the absolute error in Mv is sim- 
ply given by a = e x |My|. As shown in Bonatto & Bica 
(2011), the end result is a distribution function significantly 
smoother than a classical histogram. Indeed, the observed 
LF (Fig. 1) is approximately gaussian-shaped with a pro- 
nounced peak around Mv ~ — 7 and « 4 mags of full width 
at half maximum. 

As usual with respect to distribution functions, the 
observed (j)(M v ) can be linked to the PDMF, £(M) = 
dN/dM, by 



4>(M V ) = £(M) 



dM 



(1) 



dM v 
However, except for a few cases, the stellar mass of most GCs 
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has not been directly measured, and so, the actual shape of 
£,{M) remains uncertain. On the other hand, (,{M) can be 
inferred if one knows the initial 1 MF, the GCs age, and the 
mass-loss physics. Formally, if (, (M ) is the MF at to - and 
assuming that GCs are not formed subsequently, the MF at 
a later time t is given by 

0M o (2) 



Z(M,t) = MMo) 



dM 



where M — M(M ,t), and Mo represents the individual 
GC masses at t . 

Several works (e.g. Portegies Zwart, McMillan & Gieles 
2010; Larsen 2009) suggest that star clusters form with a 
mass distribution that follows the truncated power-law of 
Schechter (1976) 



UM ) = k.MZ 2 e-( M '' M °) 



(3) 



over the mass range [Minf ~ M SU p\, where M c is the expo- 
nential truncation mass, and k 3 is a normalisation constant. 
We note that, because of the continuous mass loss (see be- 
low), GCs that formed with a relatively low mass may no 
longer be detected, and can be considered as dissolved. Here, 
we represent this dissolution limit by the lower observable 
mass value Mi- 

The earliest cluster phase is dominated by effects re- 
lated to the impulsive removal of the parental, intra-cluster 
gas by supernovae and massive-star winds, which rapidly 
damps cluster formation and changes the ratio between 
residual gravitational potential and stellar velocity disper- 
sion, thus leading to the first bout of dynamical evolution 
(e.g. Lada, Margulis & Dearborn 1984; Goodwin 1997; Geyer 
& Burkert 2001; Baumgardt & Kroupa 2007). Another con- 
sequence of these early effects on the ICMF may be a deple- 
tion in the number of GCs towards the low-mass range (e.g. 
Kroupa & Boily 2002; Baumgardt, Kroupa & Parmentier 
2008; Parmentier et al. 2008). Clusters can be considered 
to be essentially gas free after the first t ~ 10 7 yrs, but 
they keep losing mass mainly through stellar evolution and 
externally-driven tidal effects. After several Gyrs, these pro- 
cesses may reduce the mass to a fraction of the initial value. 
Following Lamers et al. (2005), we express the mass remain- 
ing (due to stellar evolution and tidal effects) in a given GC 
(with initial mass Mo) at time t by 

l/ 7 

■yt ( Mo 



M(t) = Mo 



/',, 



1O 5 M 



(4) 



where /i se is the stellar evolution mass loss, t$ is the dis- 
solution timescale of a 10 5 Mq star cluster, and 7 is the 
power-law exponent that sets the dependence of the cluster 
disruption time on mass. We note that fi se is a function of 
time, but since we are dealing with GCs with ages between 
9-13 Gyr, we simplify this point by taking here the asymp- 
totic value /Use ~ 0.7 (Lamers et al. 2005). Useful Milky 
Way approximations for the other parameters are 7 ~ 0.7 
(Lamers, Baumgardt & Gieles 2010), and T5 « 10 4 Myr 
(Kruijssen & Cooper 2011 and references therein). How- 
ever, r 5 is expected to present a significant dispersion around 

1 In the present context, initial refers to the onset of the gas- free 
evolution. Hereafter, we refer to the beginning of this period as 

to- 



its quoted value to account for different GC environments 
and orbital conditions (Kruijssen & Cooper 2011). In this 
sense, we consider T5 to be a free parameter allowed to vary 
within a given range (Sect. 3). Thus, inverting Eq. 4 to have 
Mo = M a (M(t)) and coupling Eqs. 2 and 3, the PDMF 
can be expressed as 



£{M,t) = ksHseM"' - L {M~< + pt) 



^-(1+1/7) e - 



(MTI+fit) 1 /-! 



,(5) 



where f3 = (-yl0 5 "')/T 5 . 

Formally, Eq. 5 applies to an individual cluster having 
evolved over a time t after t . However, the underlying idea 
of our approach is to consider not the individual, but the 
collective - or average - evolution of the full sample of the 
Galactic GCs until now. This adds another component to 
our approach because, instead of a single formation event, 
the Milky Way GCs present some dispersion in age (e.g. 
Marin- Franch et al. 2009). In this sense, t should be taken 
as the mean - or at least, the representative - age (t gc ) of 
the GC population. So, hereafter we consider t — t gc . The 
constant k s can be computed by noting that the integral of 
Eq. 5 over the minimum observable GC mass, Mi, and the 
maximum remaining mass at time t gc , M s = M S u P (t gc ), 
corresponds to the present-day number of GCs (Ngc) 

N G c= / £{M,t gc )dM. (6) 

Changing variables in the integral, Eq. 6 can be reduced 
to 



Ngc 



[rln(r)] dr 



M c I ' M, 

where El is the exponential integral, 

j 7 \ 1 /t 



£l(ln(r)) 



rln(r) 



ln(Ti) = (MJ +/3t gc ) ''< I {n se M c ) , 



and 



ln(T.) = (Ml +pt gc ) lh I ifi Be M c ) . 

Summarising, the predicted PDMF (Eq. 5) can be in- 
serted into Eq. 1 to build the theoretical LF <p(Mv) and, 
with the usual relation between luminosity (Lv) and abso- 
lute magnitude L v = W -oa(m v -m®) ^ with M & = 4 g3; the 
theoretical LF is expressed as 
dM 



<t>(M v ) = 0.4t(M,t gc ) 



d log Lv 



(7) 



Subsequently, <fi(Mv) can be compared with the ob- 
served <fi(Mv) to search for the parameters that lead to the 
best match between both. Note that the present-day GC 
masses in Eq. 7 should be expressed as a function of Mv 
(through Lv), which can be done by means of a mass-to- 
light relation (MLR = M/L v ). 

2.1 The mass-to-light ratio 

According to Sect. 1 and Eq. 7, it is clear that the depen- 
dence of MLR on luminosity (or mass) is central to the task 
of building MFs from the observed LF. In this sense, we con- 
sider here the following comprehensive range of MLR shapes, 
(i) constant: M/Lv — a; (it) linearly- increasing with lu- 
minosity: M/Lv — a + bLv; (Hi) power-law: M/Lv — 
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Figure 1. The observed (light-shaded symbols) and best-fit LFs 
(solid line) produced with alternative MLRs. The respective Rrms 
are shown. A constant MLR with the mean value of M/Ly = 2 
(dotted line) taken from the literature produces a relatively poor 
fit with Rrms ~ 0.36 (top-left panel). 



a + bLy\ and (iv) logarithmic: Ai/Lv = a + b log(Ly), 
where a and b are constants. Note that for practical reasons 
(i.e., the need of an analytical expression to convert My into 
M in Eq. 7), we express the MLR as a function of Ly and 
not mass. 

Obviously, once the MLR constants a and b have been 
assigned values, the minimum and maximum present-day 
GC masses (Mi and M s ) are naturally determined from the 
respective faint and bright boundaries of the observed LF. 
And, for a given mean GC age (t gc ), a similar reasoning leads 
to the maximum GC mass (Ms Up ) at t through inversion 
of Eq. 4. 



3 SEARCHING FOR THE OPTIMUM 
PARAMETERS 

The approach summarised by Eq. 7 contains some poorly- 
known (and unknown) parameters, such as the mean cluster 
age (tgc), the dissolution timescale (£5), the truncation mass 
(M c ), and the MLR constants a and 6. Our approach consid- 
ers these five parameters as free. However, we note that some 
dependence among these parameters should be expected to 
occur, especially for very wide search ranges. For instance, 
a short t gc would result in a low t$ to hasten dissolution 
over a shorter age range. Similarly, a lower MLR would shift 



the MF turnover to a lower cluster mass, which could result 
from a longer dissolution timescale (i.e. less dynamical evo- 
lution) for a given mean cluster age. We minimise this effect 
by adopting a reasonable GC age range (see below). 

Thus, the task now is to search for values for the set 
of free parameters P = {a, b, t gc , £5, M c \ that produces 
the best match between the observed and theoretical LFs, 
4>{My) and <f)(My), respectively. This is achieved by min- 
imising the root mean squared residual (Rrms) between both 
functions. By construction, R rms is a function of the free pa- 
rameters, i.e., Rrms = Rrms(P), which we express as 



Rr 



\ 



N G c 



E 



[4>{M: 



<t>(My 



4>(M V ) 



+ HMy) 



(8) 



where the sum occurs over the nj non-zero magni- 
tude bins of the observed LF. The normalisation by the 
theoretical+observed density of GCs in each bin gives a 
higher weight to the more populated bins 2 . Finally, the 
squared sum is divided by the total number of observed 
GCs, which makes Rrms dimensionless and preserves the 
number statistics when comparing LFs built with unequal 
GC populations. 

Among the several minimisation methods available in 
the literature, the adaptive simulated annealing (ASA) is 
adequate to our purposes, because it is relatively time effi- 
cient, robust and capable to distinguish between different lo- 
cal minima (e.g. Goffe, Ferrier & Rogers 1994). The first step 
in the minimisation process is to define individual search 
ranges (and variation steps) for each free parameter. Next, 
initial values are randomly selected for all free parameters 
(i.e., the initial point P in the Rrms hyper-surface) and the 
starting value of Rrms is computed. Then ASA takes a step 
by changing the initial parameters (P — T P), and a new 
value of -Rrms is evaluated. Specifically, this implies that a 
new LF is built with the changed parameters. By defini- 
tion, any step that decreases R r ms (downhill) is accepted, 
with the process repeating from this new point. However, 
uphill steps may also be taken, with the decision made by 
the Metropolis (Metropolis et al. 1953) criterion, which has 
the advantage of enabling ASA to escape from local min- 
ima. The variation steps become smaller as the minimisa- 
tion is successful and ASA approaches the global minimum 

The output of a single ASA run is the set of optimum 
values Pmin (expected to correspond to the global minimum) 
and the respective R rm s , but with no reference to uncertain- 
ties. However, given the errors associated to the observed LF 
(Fig. 1), uncertainties to the derived parameters are strongly 
required. Thus, we run ASA several times (N rU n) allowing 
for different starting points (P ) and taking the uncertainty 
of each My bin in the observed LF (Fig. 1) into account. In 
the end, we compute the weighted mean of the parameters 
over a range of runs, using the Rrms of each run as weight 
(w — 1/Rrms)- As a compromise between statistical signifi- 
cance and computational time, we adopt here N rU n = 6000. 



2 In Poisson statistics, where the uncertainty (cr^,) of a signal <f> 
is o^ = \l <j>, our definition of Rrms turns out being equivalent 



to the usual \ 
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Free Parameters 



Derived Parameters 



rtanK firms gc 

(Gyr) 
(1) (2) (3) 



(-1) 



(5) 



M c 

(M ) 

(6) 



Mi 



M s 



(Gyr) 
(7) 



(1O 3 M ) (1O 7 M ) 
(8) (9) 



(io 4 ) 

(10) 



M °GC 



1V1 GC 



(1O 8 M ) (1O 6 M ) 
(11) (12) 



Constant: M/ Ly 
(xlO 5 ) 



Abs.Min. 0.1799 10.7 0.59 

1% Low. 0.1875 11.4 ±1.0 0.49 ±0.19 
Overall 0.2041 11.3 ± 1.0 0.50 ± 0.23 



6.0 5.8 

5.1 ±2.0 8.1 ±4.3 
5.5 ±2.7 8.7 ±5.3 



0.7 0.5 9.9 4.4 9.2 

0.7 ±0.3 0.4 ±0.2 8.4 ±1.7 3.7 ±1.4 8.3 ± 3.2 
0.8 ±0.4 0.4 ±0.3 7.5 ±2.0 3.4 ±1.5 8.4 ± 3.9 



(xlO- 1 ) 



Logarithmic: M/Ly = a + b log(Ly-) 
(xlO 6 ) 



Abs.Min. 0.1667 10.3 0.32 0.19 1.6 5.0 0.7 

1% Low. 0.1772 11.1 ±1.1 0.19 ±0.21 0.15 ±0.05 1.3 ±0.5 7.8 ± 4.5 0.8 ± 0.3 
Overall 0.1936 11.4 ±1.0 0.65 ± 0.66 0.16 ±0.05 1.5 ±0.6 6.9 ± 2.7 0.9 ± 0.3 



1.2 7.3 3.8 17 

0.7 ±0.4 5.8 ±1.2 2.9 ±1.0 14 ± 5 
0.8 ±0.5 5.3 ±1.1 3.1 ±1.0 16 ±5 



Linear: M/ Ly 
(xl0~ 5 ) (xlO 10 ) 



bL x 



Abs.Min. 0.1298 
1% Low. 0.1345 
Overall 0.1543 



12.6 0.43 0.20 1.6 

11.4 ±1.1 0.91 ±0.47 0.45 ±0.23 5.2 ± 2.8 
11.3 ±1.1 0.93 ±0.44 0.46 ±0.22 5.0 ± 2.8 



5.4 


0.6 


26 


1.9 


1.2 


21 


9.7 ±6.1 


1.4 ±0.8 


13 ±15 


1.7 ±0.6 


2.7 ±1.4 


41 ±22 


9.1 ±5.7 


1.5 ±0.8 


18 ±22 


1.7 ±0.3 


1.7 ±1.3 


42 ±20 



Power-law: M/Ly 
(xl0~ 3 ) (xlO 10 ) 



■bL° v 67 



Abs.Min. 0.1034 10.5 0.83 0.36 1.6 5.8 

1% Low. 0.1213 11.3 ±0.9 0.72 ± 0.33 0.29 ±0.14 5.1 ±2.8 8.8 ± 4.8 
Overall 0.1395 11.3 ± 0.8 0.71 ± 0.34 0.28 ±0.13 5.0 ± 2.8 9.0 ± 5.5 



1.8 5.3 1.7 3.1 44 

1.2 ±0.6 7.0 ±6.7 2.2 ±0.4 2.7 ±1.3 37 ± 17 
1.2 ±0.6 6.8 ±6.8 2.1 ±0.4 2.6 ±1.2 36 ± 17 



Table Notes. Cols. (1) and (2): Rank and corresponding R rm s value; Col. (3): mean GC age; Cols. (4) and (5): MLR constants; Col. (6) 
exponential truncation mass; Col. (7): dissolution timescale of a 10 5 Mq star cluster; Col. (8): lowest present-day (observed) GC mass; 
Col. (9): upper cluster mass at i ; Cols. (10) and (11): number and total mass of GCs at t (with Mi n f > 10 2 Mq); Col. (12): present-day 
total GC mass. 



The search ranges adopted here try to emulate the phys- 
ical conditions prevailing in the Milky Way. In this sense, the 
mean GC age runs through 9 < t gc (Gyr) < 13; the exponen- 
tial truncation mass within 10 4 ^ _M C (M ) ^ 10 11 ; and the 
dissolution timescale within 5 ^ rs(Gyr) ^ 30 (to account 
for different environments and orbital conditions). Initially, 
the MLR constants a and b were allowed to vary within 
and 10 but, after a few runs, their ranges were fine-tuned to 
more constrained ranges. 



4 RESULTS AND DISCUSSION 

The approach described in Sect. 3 was applied to the four 
MLR modes discussed in Sect. 2.1, and the numerical re- 
sults of the 6000 independent ASA runs are summarised 
in Table 1. For conciseness, we show only the parameters 
belonging to the run with the absolute minimum Rrms, to- 
gether with those averaged over all the 6000 runs, and those 
corresponding to the average of 1% of the runs having the 
lowest Rrms- The LFs corresponding to the absolute mini- 



mum Rrms are compared to the observed LF in Fig. 1. Given 
the relative proximity among the minimum and average pa- 
rameters (Table 1), the respective LFs of each MLR should 
also be similar. 

Qualitatively, the adopted MLRs produce similarly 
evolved LFs (Fig. 1), but with different individual and inte- 
grated parameters at the output. As expected, the mean age, 
in all cases, is consistent with the age of the Galaxy and the 
age spread of the Galactic GCs (Muratov & Gnedin 2010), 
ranging within 10 and 12.6 Gyr. Also, the best-fit dissolu- 
tion timescale falls in the range 5 ^ Ts(Gyr) ^ 10, a spread 
that is consistent with the different environments where the 
Galactic GCs dwell. Below we discuss the results for the 
MLR modes in a decreasing order of R rm s ■ 

Before doing so, it is useful to gather some integrated 
Galactic GC values to be used for comparison with our re- 
sults. For instance, based on the uniform ratio M/Ly = 2, 
Mackey & van den Bergh (2005) estimate that the current 
total mass of halo GCs as 2 x 1O 7 M . However, this may 
be a lower limit to the actual value of Mqq" , since halo 
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Figure 2. Present-day (empty circles) and initial (dashcd-linc) 
MFs that result from different MLRs. The PDMFs (expressed as 
dN/d log(.M)) are fitted with lognormal curves (thin solid line). 



GCs correspond to a sub-population of the Galactic GCs 
(e.g. Mackey & van den Bergh 2005) and SSP models (e.g. 
Bruzual & Chariot 2003; Anders & Fritze-von Alvensleben 
2003) usually predict the higher ratio M/Ly =2 — 4. 
Also, the initial total mass in GCs (M GC ) can be com- 
pared with estimates of the current stellar halo mass of 
M h ~ 1 x 10 9 M Q (e.g. Freeman & Bland-Hawthorn 2002). 

Constant [M/L v = a}: With R rms ~ 0.18 - 0.2, the 
best fit was obtained with M/Ly « 0.6. For such MLR, the 
lowest present-day GC mass turns out being somewhat low 
for a GC, Mi ~ 700 M ; the exponential truncation mass is 
M c ~ 6 x 10° Mq. Reflecting the low statistical significance 
of the constant MLR, the dispersion around the mean for 
most parameters is very large. When combined, the derived 
parameters imply an exceedingly-low (see above) total mass 
for the present-day GC population, M^°c < 1 x 1O 7 M . 
Assuming 100 Mq as the lower-limit for cluster formation, 
the total stellar mass stored in clusters at t would be 
Mq C ~ 4 x 10 8 M Q . Thus, the fraction of the initial stellar 
mass still bound in GCs would be as low as fbound ~ 2%. 
The higher ratio M/Ly — 2 (the average used in previ- 
ous works, e.g. Kruijssen 2008) leads to an evolved LF that 
fails to describe the observations, especially around the peak 
(My ~ — 7), with the high value of Rrms ~ 0.36. 

Logarithmic [M/Ly = a + b log(Ly)] : With R rm s 
similar to the constant MLR, the logarithmic MLR yields 
M t ~ 800 M , M c ~ 1.5 x 1O 6 M , Mg° c w ~ 1.5 x 1O 7 M , 
and M GC ~3x 10 8 M , which implies the low bound mass 
fraction of fbound ~ 5%. These values are similar to those 
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Figure 3. The best-fit MLRs are shown as a function of the clus- 
ter mass. For comparison, the curves corresponding to the abso- 
lute minimum (solid line) and the 1% runs with the lowest Rrms 
values (dot-dashed) are shown, together with the la uncertainty 
domain of the latter. 



obtained with the constant MLR and, as such, may not be 
realistic. 

Linear [M/Ly = a + b Ly] : Dropping somewhat to 
Rrms ~ 0.13 — 0.15, the linearly increasing MLR pre- 
dicts the more realistic values of Mi ~ 1.4 x 10 3 M and 
Mg°S" ~ 4 x 10 7 M . With M% c ~ 2 x 10 s M Q , the fraction 
of stellar mass still bound in GCs rises to fbound ~ 20%. 
Interestingly, the best-fits require an exponential truncation 
mass of M c ~ 10 10 M , about 4 orders of magnitude higher 
than usually adopted (e.g. McLaughlin & Fall 2008). 

Power-law [M/Ly = a + bLy] : To minimise the num- 
ber of free parameters, we applied our approach for sev- 
eral fixed exponents n, with the lowest values of R rm s 
(~ 0.10 - 0.14) obtained with n = 0.67. The relevant best- 
fit parameters are Mi ~ 1O 3 M , M c ~ 5x 1O 1O M , 
Mg°S" ~ 4 x 10 7 M , and M% c ~ 3 x 10 s M , with a bound- 
mass fraction of fbound ~ 15%. 

Both the linear and power-law models predict values 
of Mq™ that are compatible with independent estimates. 
They also imply an initial total GC mass that corresponds 
to < 30% of the current stellar halo mass. 

Reflecting the (cluster) mass-dependent nature of the 
mass-loss processes, GCs having a mass within (1 — 2) x 
10 5 M at the onset of the gas- free phase would have evolved 
to become the present-day least massive {Mi ~ 1 x 10 3 M ) 
GCs. This means that, after a Hubble time of evolution in 
the Galaxy, GCs occupying the low-mass tail of the PDMF 
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Table 2. Lognormal parameters of the present-day MF 
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<M ) 
(3) 
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(1O 5 M ) 
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(6) 


Constant 


3.5 ± 0.2 


0.80 ± 0.02 


0.970 


0.6 


1.6 


Logarithmic 


5.7 ± 0.4 


0.86 ± 0.02 


0.989 


1.0 


1.8 


Linear 


3.1 ± 0.2 


0.90 ± 0.01 


0.997 


2.7 


8.6 


Power-law 


7.8 ± 0.2 


0.91 ± 0.01 


0.997 


2.5 


3.2 
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log(AI); 


Col. (4): 
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retain only ~ 1% of their initial stellar mass. Regarding 
the total stellar mass still bound in GCs, the best-fit MLRs 
(power-law and linear) imply that ~ 85% of the mass of the 
initial GC population would have been lost to the field. 

Finally, the best-fit MLRs are shown in Fig. 3 as a func- 
tion of the cluster mass. Both the constant and logarithmic 
MLRs are somewhat shallow, with M/Lv < 2.0 over the 
full range of mass (10 2 < M(M ) < 10 7 ). On the other 
hand, the linear and power-law MLRs have M/Lv ~ 1.0 
for 10 2 ^J A4(Mq) < 10 4 , then rising monotonically for 
larger masses, reaching M/Lv ~ 2.5 for M — 10 6 Mq, and 
M/L v ~ 4 - 8 for M — 10 7 M Q . When converted to a func- 
tion of GC mass, the power-law MLR follows the relation 



M/Lv = (0.78±0.01) + 



M 



(2.17 ±0.06) x 10 5 M Q 



In summary, the best results in terms of minimum R rm s 
residuals and realistic GC parameters are obtained with 
MLRs that are relatively low for most of the GC mass range, 
but increase with luminosity (or mass), in the present case, 
expressed as a power-law (exponent 0.67) or linearly. 



ejection 4 , Eq. 4 probably yields a softer mass- loss rate than 
those previously used. We also note that, for masses larger 
than that of the turnover, the PDMFs (of the power-law and 
linear MLRs) have a slope similar (within the uncertainties) 
to the MFs of young clusters and molecular clouds observed 
in the Milky Way and other galaxies (McLaughlin & Fall 
2008, and references therein). A similar high- mass slope was 
found for a sample of Galactic and M31 GCs by McLaugh- 
lin & Pudritz (1996). Actually, this is a natural consequence 
of the large M c , which favours the formation of a higher 
number of massive clusters with respect to a low-A'fc ICMF. 
Massive clusters are naturally long-lived and, thus, the slope 
of the large-mass tail of the ICMF tends to be preserved for 
longer periods of time. 

We note that our best-fit MLR model predicts a maxi- 
mum PDMF mass (Fig. 2) that is compatible with current 
estimates for Omega Cen (NGC5139): 2.5 x 10 6 M Q (van de 
Ven et al. 2006), 4 x 10 6 M Q (Pryor & Meylan 1993), and 
5.1 x 10 6 M Q (Meylan et al. 1995). In particular, the latter 
estimate implies the ratio M/Lv ~ 4.1, again consistent 
with the power-law MLR value for a GC with mass in the 
range (2 - 5) x 10 6 M Q (Fig. 3). The power-law MLR pre- 
dicts the mean present-day GC mass M ~ 2.5 x 1O 5 M0, 
from which follows the scaling relation M ~ 3Mto- 

We also show in Fig. 2 the ICMFs reconstructed with 
the best-fit parameters. Interestingly, while the ICMFs pre- 
dicted by the constant and logarithmic MLRs show the typ- 
ical break at M c ~ 10 5 Mq of the Schechter profile, those 
of the linear and power-law MLRs do not. Both are charac- 
terised by the large value of M c ~ 10 M© and, thus, they 
end up resembling simply as a scale-free, power-law of slope 
— 1 (or —2, when expressed as dN/dM) over the very wide 
mass range 10 2 — 10 10 Mq. The latter is a typical feature of 
young cluster MFs observed in nearby galaxies (e.g. Zhang 
& Fall 1999; Portegies Zwart, McMillan & Gieles 2010). 



4.1 The present-day MF 

One of the main results of our approach is to provide access 
to the PDMF through Eq. 5 and the best-fit parameters of 
each MLR mode. For practical reasons, they are shown in 
terms oidN/d log(A4) in Fig. 2. They all clearly resemble a 
lognormal mass distribution 3 , which is confirmed by the cor- 
responding fit - the parameters are given in Table 2. Again, 
the lognormal character applies especially to the MFs pro- 
duced by the power-law and linear MLRs, followed by the 
logarithmic and constant modes. 

Overall, the PDMF mass turnover predicted by our 
MLRs occur within Mto ~ (3 - 8) x 10 4 M (the upper 
bound corresponding to the power- law MLR), a range some- 
what lower than the Mto ~ (1 — 2) x 10 J M© derived in pre- 
vious works (e.g. McLaughlin & Fall 2008). Probably, the 
difference occurs because for most of the GC mass range, 
our MLR models are lower than the M/Lv — 2 usually 
employed in previous studies, which naturally shifts Mto to 
lower values, implying longer dissolution timescales. In addi- 
tion, by not including the effects of shocks with giant molec- 
ular clouds and spiral arms, and two-body evaporation and 



5 SUMMARY AND CONCLUSIONS 

We present an approach to recover both the initial and 
present-day mass functions of the Galactic globular clusters, 
having as constraint the observed luminosity distribution. 
This is achieved by taking a few mass-loss processes into 
account and assuming a mass-to-light ratio that is either 
constant or increases with luminosity linearly, as a power 
law or logarithmically. 

It starts with a Schechter-like ICMF, where the trun- 
cation (Mc) mass is a free parameter. Then, using a set 
of parameters that usually apply to the dynamical dissolu- 
tion of Galactic star clusters, among these the mean GC 
age (tgc), the dissolution timescale of a 10 Mq cluster (ts), 
and an MLR (parametrised by the constants a and b), we 
compute the MF shape after evolving for a time t gc after 
the onset of the gas-free phase. The evolved MF is subse- 
quently converted into the LF. Finally, we search for the op- 
timum parameters (tgc, is, Ale, a, b) that minimise the resid- 
uals (Rrms) between the evolved and observed LFs. 

The best results, both in terms of minimum values of 
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4 Although molecular clouds and spiral arms may not be relevant 
to most GCs, evaporation does play a role (e.g. McLaughlin & Fall 
2008). 
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Rrms and realistic parameters for the Galactic GCs, cor- 
respond to the power-law and linear MLRs, respectively. 
Specifically, because both MLRs increase with luminosity 
(or mass), their present-day MFs imply a total mass in GCs 
of the order of Mq°c ~ 4 x 10 7 M , which corresponds to a 
fraction of ~ 15% of the initial stellar mass in GCs. On the 
other hand, the constant and logarithmic MLRs are exceed- 
ingly shallow for the full GC mass range, yielding unrealis- 
tically low values of M™q . 

Because of the large truncation mass (M c ~ 10 10 Mq), 
the Schechter-like ICMFs (of the linear and power-law 
MLRs) end up resembling the scale-free, power-law MFs 
of young clusters and molecular clouds in the Milky Way 
and other galaxies. In addition, their PDMFs follow closely 
a lognormal distribution with a turnover mass at Mto ~ 
7 x 10 4 Mq that, for larger masses, behave as a power-law 
with the same slope as the ICMFs. The low values of Mto, 
compared to previous works, follow from the relatively low 
MLRs (for most of the GC mass range) dealt with here. 

Summarising, it is clear that the power-law is not the 
unique ICMF solution, the MLR, which serves to fine-tune 
the details, must increase with cluster mass (or luminos- 
ity), and that mass loss is the main factor responsible for 
shaping the ICMF into its present, lognormal form. Finally, 
our results suggest a common origin - in terms of physical 
processes - for GCs and young clusters. 
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